Three-party pure quantum states are determined by two two-party reduced states 
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We can uniquely calculate almost all entangled state vectors of tripartite systems ABC if we 
know the reduced states of any two bipartite subsystems, e.g., of AB and of BC. We construct the 
explicit solution. 
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Generic multiparty composite quantum states estab- 
lish complex multiparty correlations. In the particular 
case of pure composite states, however, recent evidences 
have shown that higher order correlations follow from 
lower order ones 0, 0, Such quantum features 3 
are of central interest in the modern field of quantum in- 
formation H as well as in the more traditional field of 
many-body physics 0. 

Generic three-party pure quantum states have shown 
to be uniquely determined by their two-party reduced 
states 111 |2j. Consider, e.g., a composite pure state 
Pabc = \i>ABc){i>ABc\ of three parties A,B,C of di- 
mensions dA,d,B,dc, respectively. Let pab, Pbc, Pac 
denote the two-party reduced states. In case of three 
qubits, these three reduced states will uniquely deter- 
mine the composite state pabc in almost all cases pj. 
For higher dimensions, satisfying the 'triangle-inequality' 
dA > ds + dc — 1, an alternative theorem holds: the two 
reduced states pab and pac are already sufficient to cal- 
culate the state pabc of the whole system . Note that 
in both cases one calculates pabc without assuming that 
Pabc is pure. It comes out from the reduced states. If 
one assumes it then a stronger statement holds. As we 
shall prove in the present work, almost all pure compos- 
ite states \ipABc) can uniquely be calculated from the 
knowledge of any two of the two-party reduced states if 
one knows already that pabc is pure. This result holds 
in any finite dimensions. We present explicit equations 
for \4>abc)- 

For concreteness, let us prove how a generic \tpABc) is 
determined by pab and psc- Obviously, the latters de- 
termine the three single-party reduced states pa,Pb,Pc 
as well. One shall diagonalize them, e.g.: 

PA=J2pA\i)(il Pa>0- (1) 

i 

Similarly, \j) and |fe) stand for the eigenvectors with non- 
zero eigenvalues p B ,Pc of ps and pc respectively. Since 
Pabc is pure, the reduced state pa shares its eigenvalues 
p\ with pbc- 

PBC = J2PA\i;BC)(i;BC\, (2) 

i 

where \i\BC) are the orthogonal eigenvectors of pbc 



with non-zero eigenvalues. Similarly, we introduce the 
orthogonal decomposition of pab as well, with non-zero 
eigenvalues p k c and eigenvectors \k;AB). We may omit 
decomposition of pac'- it is not required by the present 
proof. From the spectral decompositions (1,2) we can re- 
construct the Schmidt-decomposition of all three-party 
pure states compatible with pa and pbc'- 

\tpABc; a) = ^2 ex P( ia 0\/p^l*) ® K; BC), (3) 

i 

where a = is the set of phases to be specified later. 
From the spectral decompositions of pab and pc, we 
have another family of all pure states compatible with 
Pab and p c : 

Habc-,7) =5>xp(i 7fe ) v /z£|fc;AB)®|fc). (4) 
fe 

Since the true \ipABc) is compatible with both pab and 
Pbc (and thus with pa, pc) therefore at least one solution 
exists for the a^s and 7fe's such that 

\ipABC;a) = \ipABC;j)- (5) 

We are going to prove that this solution is unique hence 
the state (5), derived from pab and pbc, will be the true 
\iPabc)- 

First we cast the vectorial equation (5) into equations 
for amplitudes. Let us calculate the following coefficients: 

A) k = (jk\i;BC), C^ = {ij\k;AB), (6) 

They are non- vanishing for a generic state \ipABc)- In 
fact, the eigenvectors (with non-zero eigenvalues) of a 
composite state are superpositions of the direct-products 
formed by the eigenvectors (with non-zero eigenvalues) of 
the respective subsystem reduced states 0- In our case, 
we use the following expansions: 

\i;BC)=J^A) k \jk), \k;AB) = J^C*\ij). (7) 

jk ij 

Substituting them into eqs.(3) and (4), considering or- 
thogonality of the product states \ijk), we expand eq.(5) 
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into the following set of compatibility equations between 
a and 7: 

exp{ia t )^f~A) k = exp(ii k ) iJifeCtj (8) 

for all i,j, k. Multiplying the l.h.s by the complex conju- 
gate of the r.h.s. and the r.h.s. by the complex conjugate 

of the l.h.s. will cancel the factors J PaPc > yielding: 

exp[i(a, - 7fc)MjfcCy- = exp[-i(a, - lk )]X ]k C^. (9) 
Finally we obtain the following simple equations: 

a i -7fe = arg^^ fe C^. (10) 

3 

for all i and k. The solution cti,"f k is then trivial and 
unique upto an (irrelevant) constant phase shift — > 
a i + X,lk ~~ * 7fc + X- The constant x contributes to an 
irrelevant phase factor exp(i^) in front of the pure state 
(5). 

Ref.]3| considers the generic pure state \tp) of a large 
number of identical parties of dimension d each. The au- 
thors derived the upper bound ajj = 2/3 on the fraction 
of parties whose reduced states enable one to reconstruct 
\ip). The lower bound = 1/2 was obtained for large 
d. It should be observed that, for the lower bound, the 
authors assume that one restricts the reconstruction for 
pure states. For such conditions, our alternative theorem 
will sharpen the upper bound ajj = 2/3. Let us group 
the parties into three subsystems A,B,C where, e.g., A 
is a single ci-state system while B and C share on the 
rest equally or almost equally. According to our theo- 
rem, pab and pac determine a generic pure state of 
the whole system. This yields ajj = 1/2 asymptotically. 
Observe the coincidence with the lower bound oll = 1/2 
Q. Accordingly, there must be an (almost) one-to-one 
mapping between the space of pure states of the whole 
system and (a certain region in) the space of the reduced 
states of all fractions ~ 1/2 of the whole. 

We summarize the steps reconstructing a state vec- 
tor \4>abc) from two density matrices pab and psc- 
First we calculate pa, Pb, pc- Then we diagonal- 
ize pab, Pbc, PA, Pb, pc, and calculate the coefficients 
A l j k ,Cfj (6). The wanted pure state \i/jabc) takes the 
form (3) with 

a i = axg s ^J! jk Cl , (11) 

3 



where k is set to any fixed value. Recall that the Aj k s 
and C^-'s are not independent at all. The above particu- 
lar expression of could well be replaced by a variety of 
equivalent expressions of them. One can, e.g., take any 
fixed value for j instead of the summation over j which 
we took for representation invariance. To display explicit 
invariance of the full solution we need further investiga- 
tions on the underlying geometric structure. 

Finally we mention a possible extension of the method 
for spatial tomography. Assume that we have to recon- 
struct a spatial wave function i\){xyz) from planar pro- 
jections. Let us define the density matrices in the XY- 
and FZ-planes, e.g.: 



PxY{xy;x'y') = j i/}(xyz)ip(x'y' z)dz, (12) 



and a similar equation for pyz- If our theorem re- 
mains valid for infinite dimensions as well then the re- 
duced states pxy, Pyz would determine the original spa- 
tial wave function ?p(xyz). The choice of concrete equa- 
tions of reconstruction would then require special care 
against numeric instabilities. 
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